in 
< 

(N 



The hole spectral function and the relationship between overlap 
functions, natural orbitals and the one-body density matrix in 

nuclei IjThis work is partially supported by the Bulgarian National 
Science Foundation under the Contracts Nr. $-32 and $-406 and by the 
Royal Society and the Bulgarian Academy of Sciences. 

A.N. Antonov, M.V. Stoitsov, M.K. Gaidarov, S.S. Dimitrova 

Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Sofia 

1784, Bulgaria 

and 



> 

00 ! P-E. Hodgson 

^ ! Nuclear Physics Laboratory, Department of Physics, University of Oxford, Oxford 

00 
O 

vn 
0^ 



0X1-3RH, U.K. 

A method to calculate the hole spectral function in the discrete part of the spectrum is suggested 
within the natural orbital representation of the one-body density matrix of ^-nucleon system using its 
relationship with the overlap functions of the eigenstates in the {A — l)-nucleon system. 



, The cross-section of direct nucleon removal processes is determined by the spectral function which contains 

■ the information on the nuclear structure and is interpreted as the probability for the removal of a nucleon 

with given momentum and energy from the target nucleus with A nucleons (e.g., [1-25]). In particular, in 
the plane-wave impulse approximation the cross-section for the direct knock-out process is proportional 
to the diagonal element ^(k, k, E) = S'(k, E) of the hole spectral function (or matrix) in the momentum 
representation: 

5(k, k'; E) - (*o a+{k')5{E + H- £;^)a(k) , (1) 



where j^'o) is the ground state wave function of the target nucleus with A nucleons, a+(k') and a(k) 
are creation and annihilation operators for a nucleon with momentum k' and k, respectively, H is the 
Hamiltonian of the system with {A — l)-nuclcons and E'^ is the ground state energy of the target nucleus. 
If the latter has a total spin and parity J"^ = 0+, then introducing a complete set of eigenstates of H for 
the system of {A — l)-nucleons j^*/) (where the state j^*/) is characterized by the energy Ef with both 
discrete and continuous values and by other discrete and continuous quantum numbers) the hole spectral 
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function can be written in the form: 



S{Kk',E) - ^{^o\a+{k')\^f){^f\aik)\^o)SiE + Ef-E[ 



(2) 

^ $;.(k')$/(k)<5(£; + i?/-£;0), (3) 



where 



$^.(k) = (vI/^|a(k)|vI/o) (4) 

is the overlap function in the momentum representation [26-28] . 

The methods used to calculate the spectral function are reviewed, e.g. in ^ |l6|. The use of 
the independent-particle shell model (when the overlap function (4) is equal to the single-particle wave 
function of the occupied state) cannot explain the fragmentation or spreading of the hole strength. This 
is because, due to the residual interaction, the hole state in the target nucleus is not an eigenstate of the 
{A — l)-nucleon system and its strength is distributed over several states of the final system. 

In this work we suggest a method to calculate the hole spectral functions using essentially the one- 
body density matrix (OBDM) of A-nucleon system in the natural orbital representation |^ and its 
relationships with the natural orbitals (which diagonalize the OBDM) and the overlap functions (4). An 
expansion of the latter in the basis of the natural orbitals is used. The following two reasons can justify 
the use of the method: 

1) Recently the diagonalization of the realistic one-body density matrix of the correlated nuclear 
ground state obtained by various correlation methods [|6), such as the Jastrow method [30-32], as well 
as the generator coordinate method [|l^, ^ Q and the coherent density fluctuation model [|l^, |3j, ^ 
gave reliable information on the natural orbitals and occupation numbers in nuclei. These quantities 
correspond to the realistic behaviour of nuclear characteristics which are sensitive to the short-range 
nucleon-nucleon correlations, such as the nucleon and cluster momentum distributions, the mean kinetic 
and removal energies, radii and others. The natural orbitals in nuclei, as well as those in other fermion 
systems, such as '^He liquid drops [ p6[ , are strongly localized and quite different from the overlap functions 
and from the mean- field type orbitals [^, ^ Thus, it is of importance to apply the natural 
orbitals corresponding to realistic OBDM obtained in correlation theoretical methods to calculate the 
hole spectral function 5(k, k',i?). 

2) The basic quantity which is necessary to calculate the spectral function (3) is the overlap function 
(4). We show in this paper that the hole spectral function in the discrete part of the spectrum can be 
calculated by using the general relationship |]37| which connects the asymptotic behaviour of the one-body 
density matrix with the overlap functions of the {A — l)-particle system eigenstates. This relationship is 
of general importance because it enables one to obtain quantities connected with the bound eigenstates 
of the (A — l)-particle system (such as overlap functions, spectroscopic factors and separation energies) 
by means of the exact OBDM (or by a realistic one obtained in a given correlation method) of the ground 
state of the A-particle system. In this way, the hole spectral function in the discrete part of the spectrum 
can be, in principle, calculated on the basis of the OBDM of the A-particle system. 

Now we introduce the necessary quantities which are used in the theoretical method to calculate the 
hole spectral function. 
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The one-body density matrix of the ground state 1^*0) of the A-nucleon system has the form 

= (^o |a+(a;)a(x')| *o), (5) 

where x = {rcrr} labels spatial, spin and isospin coordinates and {x)^a{x') are the creation and 
annihilation operators. 

The natural orbitals (NO) ipa{x) are defined |^ as the complete orthonormal set of single-particle 
wave functions which diagonalize the OBDM: 

p(x,x')=^iVa¥':(a:)v'a(a:')- (6) 

a 

The eigenvalues Na (0 < Na < 1, ^^-^a = A) are the natural occupation numbers. We note that the 

a 

sum (6) is over the discrete states determined by the finite-range NO (/3a (k). 
The OBDM (5) can be presented also in the form: 

p{x,x') = ^ 'i>}{x)^f{x'), (7) 

where ^f{x) — f \ 5'o) is the overlap function in the coordinate representation. 

The overlap functions can be expanded in the basis of the natural orbitals (e.g., in momentum space): 

$/(k)=^(^a|$/)^,(k) (8) 

a 

The hole spectral function is then given by the expression: 



aM 7 



(9) 



a,b 

where 



SabiE) EE ^ {<!>f\^a){^b\'i>f)S{E + Ef- ). (10) 

The quantity (for which different notations exist, e.g. jlj, [Ts]): 

OaJ^S'J]^{iPa\<S>f) (11) 

from (8) and (9) is the amplitude of the contribution of the orbital a to the overlap function for the 
eigenstate 1^"/). We mention that the quantity (11) determines both the spectroscopic factor of the state 

sf'' EE ($^1$^) = E ^"^/ - E K^-l'J'/)!' (12) 

a a a 

and the occupation probability of the orbital a: 

Na^i Kjf - i Saj - K^a|$/)r • (13) 

In general, for a given orbital a, only a limited subset of states / of the residual nucleus contribute to 
the sums (10) and (13). 
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The function Sab{E) (given by Eq.(lO) and often called also "spectral function") can be rewritten 
[|j introducing the different states j^*/) of the residual nucleus: i) the bound states l^* with energy 
El, and degeneracy quantum number a, and ii) the continuum states |^_e^,c} with energy Ef and the 
channel index c which specifies the channel where there is an incoming wave (all other channels contain 
only outgoing waves), as well as all degeneracies like spin projections etc. Then Eq.(lO) becomes: 

SabiE) - Y.('^-"\'P-)('Pb\<i>uc.)S{E + E,~El) 

+ Y.('^E,=E0^^EM{^b\^Ef=E0^-E,c)9iE"A " " E) (14) 

c 

EE St,'^-{E) + S:,^-{E), 

where and ^Ef,c are the overlap functions associated with the bound and continuum eigenstates of 
the residual nucleus and -E^'^T'i is the threshold for particle decay of this nucleus. If the latter is a nucleon 
threshold, then E'^J^\ = E'^_2, where i?^_2 is the ground-state energy of the nucleus with A — 2 nucleons 
[H. The hole-spectral function (14) contains two parts: i) the spectral function in the discrete part of 
the spectrum S^j^'''{E), and ii) the spectral function in the continuum of the hole spectrum S2'|^^{E) with 
E <E%- 

As can be seen from Eqs.(9) and (14) the hole spectral function is essentially connected with the 
natural orbitals {(fa} and the overlap functions and their relationship with the OBDM. Now we shall 
outline briefly this relationship. 

In the case of spherical symmetry the OBDM can be written in the form: 

pix,x') = ^p(«'^)(r,r')5]l^*„(f^,a)yy™(17',a'), (15) 

qlj m 



where the radial part of the OBDM is: 



{r,r')^^ <S>'^/'\r)'S>'^/'\r'). (16) 



In the above equations $^'''"'''(r) is the radial part of the overlap function, Yijmi^,<j) is the spin-angular 
function, q denotes the nature (proton and neutron) of the overlap function and /, j are angular and total 
momentum quantum numbers. 

It is known that the overlap functions associated with the bound states of the (A—l)- and {A + 1)- 
nucleon systems are eigenstates of a single-particle Schrodinger equation in which the mass operator plays 
the role of a potential. Due to the finite range of the mass operator, the asymptotic behaviour of the 
radial part of the neutron overlap functions for bound states i> (labeled by v—O, 1, ... with increasing 
energy) of the {A — l)-nucleon system is given by ^ ^ : 

$to)(r) C('?'^')exp(-fc('?'J'V)/r, (17) 

where 

fc(''^) = i[2m,(i?W)-i?0)]V2. (18) 

For protons some mathematical complications arise due to an additional long-range part originating from 
the Coulomb interaction p^ , though everything from the neutron case remains valid. It is assumed in 
pq that Eq.(17) is also valid for the overlap functions corresponding to the {A — 1) continuum. 
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The asymptotic form of the overlap functions (Eqs.(17) and (18)) determines the asymptotic behaviour 
of the radial part of the OBDM p^''^'' (r, r') It is shown in |^ that at large values of r' = a ^ oo 
one can derive the lowest bound state overlap function by means of the radial part of the OBDM: 



as well as the separation energy 
and the spectroscopic factor 



= h^i'"'''^ /2mg (20) 



^ (21) 



The normalization coefficient Cq'''"''' can be obtained from the asymptotic form of the diagonal part of 



the radial OBDM: 

2 



p^«'^^(a,a) d/''> exp{^2k^"'>a)/a\ (22) 

As shown in the overlap functions for all bound states of the {A — l)-nucleon system can be con- 
structed from the OBDM repeating the above procedure. For instance, the overlap function for the next 
bound state is: 

In the case of the continuum contributions to the OBDM one can calculate the particular sum over 
the scattering channels c: ''^^[^'f'^\r, E)Cc{E)], but not the overlap function for each channel [3^ . 

c 

The method for calculating of the hole spectral function in the discrete part of the spectrum (for 
non-degenerate states v): 

S^-^ {k,k', E)^J2 <(^>bik) J2{^^\^a){^b\<^.)SiE + E,- El) (24) 
from a given theoretical correlation method consists in the following procedure: 

1) By diagonalizing the one-body density matrix of the A-particle system ground state one obtains the 
natural orbitals {(pa(k)} (e.g., as in [30-35,16]); 2) The bound-state overlap functions $1/ and separation 
energies are calculated on the basis of the one-body density matrix following the approximate method 
described above (Eqs.(19)-(23)). 3) The amplitudes of the contribution of the natural orbital a to the 
overlap function {Lpa\^u) are calculated and the results substituted in Eq.(24). The overlap functions of 
the discrete states indeed allow the corresponding part of the spectral function to be calculated directly 
from Eq.(3). However, the third step of the procedure also makes it possible to calculate simultaneously 
on the same footing the amplitudes (11) and some particular terms of the sums in Eqs.(12), (13) and (24). 
The spectroscopic amplitudes (11) can then be calculated with shell model single-particle wave functions 
and compared with those obtained from natural orbitals. It can be seen from Eqs.(l) and (5) that the 
energy integral of the hole spectral function (1) defines the one-body density matrix in the momentum 
representation 

dES{k,'k!,E)^ p{k,\i!), (25) 

-00 

where Ep is a negative quantity whose absolute value is equal to the separation energy of the A-nucleon 
system [p8[. 
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We emphasize that the method described above enables one to obtain the hole spectral function in 
the discrete part of the spectrum (i.e. the integrand of the left-hand side of Eq.(25) in the energy interval 
between E^—E^_2 and Ep) on the basis of the one-body density matrix calculated in a given correlation 
method. The knowledge of p(k, k') can give some information on the remaining part of the integrand in 
the left-hand side of Eq.(25), namely the hole spectral function in the continuum part of the spectrum 
in the energy interval between — oo and E*^ — -B^_2- 

In this paper we suggest a new theoretical method to obtain the hole spectral function in the discrete 
part of the spectrum. The method is based on the natural orbital representation in nuclear theory 
and uses essentially both the natural orbitals and overlap functions as well as their relationship with 
the OBDM. Thus the theoretical point of the method consists in the possibility of using the OBDM 
which is related to the properties of the A-nucleon system to calculate the hole spectral function which 
determines the cross-section of the nucleon removal processes and gives information on the structure of 
the {A — l)-nucleon system. The applications of the method can serve also as a test of the predictions of 
the correlated methods concerning the OBDM of the correlated ground state of the A-nucleon system. 

Our program to apply the suggested method includes two stages: i) calculations of overlap functions 
on the basis of realistic OBDM from a given correlation method and studies of their properties, and ii) 
calculations of hole spectral functions in the discrete part of the spectrum and comparison with available 
experimental data. The results from the fulfilment of this program will be given elsewhere. 
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